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Abstract 



Deep inelastic scattering data on the F 2 structure function provided by the BCDMS, 
SLAC and NMC collaborations are analyzed in the non-singlet approximation with the an- 
alytic and "frozen" modifications of the strong coupling constant featuring no unphysical 
singularity (the Landau pole). Improvement of agreement between theory and experi- 
ment, with respect to the case of the standard perturbative definition of a s considered 
recently, is observed and the behavior of the higher twist terms in the next-to-next-to- 
leading-order is found to be confirming earlier studies on the subject. 



PACS : 12.38 Aw, Bx, Qk 

in ; 

Keywords: Deep inelastic scattering; Nucleon structure functions; QCD coupling con- 
stant; NNLO level; 1/Q 2 power corrections. 

00 
O 
O 



Introduction 



At present an accuracy of the data on the deep inelastic scattering (DIS) structure func- 
tions (SFs) makes it possible to study separately the Q 2 -dependence of logarithmic QCD- 
inspired corrections and those of power-like (non-perturbative) nature (see for instance [1] 
and references therein). 

This paper closely follows those devoted to the similar studies performed in [HE]. We 
analyze DIS SF F 2 (x,Q 2 ) with SLAC, NMC and BCDMS experimental data involved [4]- 
[9| up to the next-to-next-to-leading-order (NNLO) of massless perturbative QCD. This 
latter level of accuracy has become possible to ensure given the results obtained for both 
the o$.{Q 2 ) corrections to the splitting functions (the anomalous dimensions of Wilson 
operators) [10] and the corresponding expressions of the complete three-loop coefficient 
functions for the structure functions F2 and Fl [TT] . 

As in our previous papers the function F 2 (x, Q 2 ) is represented as a sum of the leading 
twist F% ® CD (x, Q 2 ) and the twist four terms: 

F 2 (x,Q 2 ) = Ff CD {x,Q 2 )(l + } ^ . (1) 

While analyzing experimental data various corrections must be taken into account. Here 
the nuclear effects, target mass corrections, heavy quark threshold corrections and higher 
twist (HT) terms are considered. For more details we refer to [3j 112]. 



1 The BFP dataset, which was included in the analyses carried out in [3J[3], is excluded in this work since its 
influence here is negligible. 



In the present work we study the effects of the popular infrared modifications of the 
strong coupling constant, such as the "frozen" |13| and the analytic [H] ones, and consider 
the change in the magnitudes of the twist four terms. These coupling constants are free 
of unphysical singularity in the infrared region of the scale ~ A. 

As is known there are at least two ways of carrying out the QCD analysis over DIS data: 
the first one (see e.g. |15[ I16j) deals with the Dokshitzer-Gribov-Lipatov-Altarelli-Parisi 
(DGLAP) integro-differential equations |17j and let the data be examined directly, whereas 
the second one involves the SF Mellin moments and permits performing an analysis in 
analytic form as opposed to the former option. In this work we adopt the analysis carried 
out over the moments of SF F^{x, Q 2 ) defined as follows 

MpQCD/twistl/...^ = f 1 x n~2 F pQCD/t W ist2/ ... ^ q2) ^ (2) 

J 

and then reconstruct SF for each Q 2 by using Jacobi polynomial expansion method |18j - 
[20J (for further details see [31 [12] and section [4]). 

The paper is organized as follows. In Sec. 2, we review theoretical aspects of the 
analysis. Sec. 3 contains infrared modifications of the strong coupling constant and the 
procedures to implement them into the representations for the SF Mellin moments. In 
Sec. 4, the SF fit method to follow is briefly described. The results obtained in two 
different schemes — the variable and fixed flavor number schemes — are given in Sec. 5. 
followed by the conclusions. In App. A, the results of the Fractional Analytic Perturbation 
Theory (FAPT) are adapted for the SF Mellin moments. App. B shows how an additional 
contribution to the structure function comes from the "analytization" of the coupling 
constant. 



2 Theoretical aspects of the analysis 

Here we describe the theoretical background of our analysis presented for the leading 
twist part. The coefficient h&(x) of the twist-four correction is considered to be Q 2 - 
independent and its x shape is determined from the fits to the DIS structure functions. 
For a bit detailed account see also [3]. 

The leading twist DIS SF can be represented as a sum of two terms: F^ wist2 '(x , Q 2 ) = 
FJ? s (x,Q 2 ) + F^(x,Q 2 ) , the nonsinglet (NS) and singlet (S) parts. At this point let's 
introduce parton distribution functions (PDFs), the gluon distribution function fc(x, Q 2 ) 
and the singlet and nonsinglet quark distribution functions fs(x,Q 2 ) and iNs(x,Q 2 )a 

f 

f 5 (x,Q 2 ) = ^f q (x,Q 2 ) = V(x,Q 2 )+S(x,Q 2 ), 
<? 

iNs(x,Q 2 ) = u v (x,Q 2 ) - d v (x,Q 2 ) , 

where / is the number of quark flavors (up, down, strange,. . .), V(x,Q 2 ) = u v (x,Q 2 ) + 
d v (x, Q 2 ) is the distribution of valence quarks and S(x, Q 2 ) is a sum of sea parton distri- 
butions set equal to each other. 

As is known, at large values of x (x > 0.25) the contribution of sea parton distribu- 
tions S(x,Q 2 ) is negligible and, therefore, fs(x,Q 2 ) and ¥ns(x,Q 2 ) (and correspondingly 
F^ s (x,Q 2 ) and F 2 5 (x,Q 2 )) appear to have similar (^-dependence. 

In the present paper we restrict analysis to the large x region. Consequently, the 
analysis is dubbed the "nonsinglet" one (simply signaling the absence of gluons) but 
actually the data on the entire SF F2(x,Q 2 ) will be considered. 



2 Unlike the standard case, here PDFs are multiplied by x. 



There is a direct relation between SF moments ([2]) and those of PDFs 



f N s(n,Q 2 ) = [ dxx n - 2 f NS (x,Q 2 ). 
Jo 

For example, in the nonsinglet case it is found to be [21]: 

M n S iQ 2 ) = RNs(f) ■ C^ st2 (n, a s (Q 2 )) ■ i NS (n, Q 2 ) , (3) 

with a s (Q 2 ) = a s (Q 2 ) / (Air) and C^ st2 (n, a s (Q 2 )) being the Wilson coefficient functions. 
The constant Rns(I) depends on weak and electromagnetic charges and is fixed to be 
one sixth for / = 4 |21j . 

Considerations of such issues as the PDF normalization, target mass (TMC) and higher 
twist corrections (HTCs), as well as nuclear effects, remain essentially the same as in our 
previous work [3] so we refer to it for further details, though quoting some salient points. 

A starting point of the evolution is taken at relatively large values Qq. There is a 
number of reasons behind that choice, e.g., fewer heavy quark thresholds have to be 
crossed to reach a normalization point, a perturbative approach must be applicable at the 
value of Qq. Besides, impact of higher order corrections derived from PDF normalization 
conditions is the more negligible the higher normalization point is. 

The moments f]\fs(n,Q 2 ) at some Qq is a theoretical input to the analysis which is 
fixed as follows. In the fits of data with the cut x > 0.25 imposed only the nonsinglet 
parton density is worked with and the following parametrization at the normalization 
point is used (see, for example, [2"2"j [23]): 

fjVs(n,Qo) = / dxx n ~ 2 f N s{x,Ql), 
Jo 

W*,Qo) = A NS (Q 2 )(l-x) b »M(l + d NS (Q 2 )x), (4) 
where A/vs(Qo)> t>Ns(Qo) an d <^vs(Qo) are some coefficients E|. 



2.1 Strong coupling constant 

The strong coupling constant is determined from the renormalization group equation to 
an accuracy of 0(1O -5 ) (which is enough for our purposes, also we checked that for higher 
precision the results get no much better). At NLO level the latter is given by 
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where hereinafter the strong coupling constant at the leading order (LO), next-to-leading 
order (NLO) and NNLO levels of approximation will be denoted by ai^(Q 2 ) with i = 0, 1 
and 2, respectively. Furthermore, we will designate ai (Q 2 ) as a s (Q 2 ) to simplify notation, 
as it is the NNLO approximation that is of special interest in the present analysis. 
At NNLO level the strong coupling constant is derived from the following equation: 
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3 Here we do not consider the term ~ a; a Ns(Qo) m normalization of fjvs(a;, Qq), because of the cut x > 0.25. 



The expression for / looks like 
I{a s {Q 2 )) = 



2 b 1 + 2b 2 a s (Q' 
arctan 
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bi + 2fe 2 a s (Q 2 ) -V^A 
&i + 26 2 a s (Q 2 ) + V /Z A 



for / = 3, 4, 5; A > 0, 



for / = 6; A < 0, 



where A = Ab 2 — b\ and b{ = /3j//3o are read off from the QCD /?- function: 



/3 a s - f3ia s - f3 2 a s + . . . . 



(8) 



2.2 Q 2 -dependence of SF moments 

The coefficient functions C 1 ^fg st2 (n,a s (Q 2 )) are further expressed through the functions 
B$ s (n) (i = 1,2), which are known exactly [331 EE] @ 

Cftg^n, *s{Q 2 )) = 1 + a a 5$7 (n) + a 2 s B {2) s (n) + O(a^) . (9) 

The (5 2 -evolution of the PDF moments can be calculated within a framework of per- 
turbative QCD (see e.g. [211 [27]) to yield: 
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h NS {n,Q 2 ) = l + a s {Q 2 )Z^ s {n) + a 2 s {Q 2 )Z^ s {n) + 0{ai{Q 2 ) 



(11) 



and 128 
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(12) 



Here 7jvs( n ) are the factors before a s in the expansion with respect to the latter of the 
anomalous dimensions "fNs( n 7 a s) (taken in the exact form from |10j). 



3 Infrared modifications of the strong coupling 
constant 

Here, we investigate the potential of modifying the strong-coupling constant in the infrared 
region with the purpose of illuminating the problem related to the Landau singularity in 
QCD. 

Specifically, we consider two modifications, which effectively increase the argument of 
the strong-coupling constant at small Q 2 values, in accordance with [29J. 

In the first case, which is more phenomenological, we introduce freezing of the strong- 
coupling constant by changing its argument as Q 2 — > Q 2 + M 2 , where M p is the rho-meson 
mass [13] (see also [301 124j and references therein). Thus, in the formulas of the previous 
section the following replacement is to be done: 

a i i\Q 2 )^af^{Q 2 ) = af{Q 2 + M 2 ) (i = 0,l,2). (13) 

4 For the odd n values, the coefficients B$ s (n) and Z$ s {ri) (i — 1,2) can be obtained by using the analytic 
continuation [2171 [231 12"6] . 



A second possibility is based on the idea proposed by Shirkov and Solovtsov [HI [31] 
(see also recent reviews [32] and references therein) regarding the analyticity of the strong 
coupling constant in the complex Q 2 -plane in the form of the Kallen-Lehmann spectral 
representation. This approach leads effectively to additional power (^-dependence for the 
DIS structure functions (see Eq. (|B2p ). 



3.1 Naive case 



This modification is in a sense quite similar to the freezing procedure given above (|13|) . 
i.e. the one- and two-loop coupling constants a^\Q 2 ) and a~P {Q 2 ) appearing in the 
formulas of Subsec. 12.11 are to be replaced as follows: 
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where the ellipsis stands for cut terms which give negligible contributions in our analysis. 

At the one-loop level, the expression for the analytic coupling constant (|14p is very 
simple. However, at higher-loop levels, it has a rather cumbersome structure (see a recent 
paper [33] and discussions therein). Therefore, it seems to be simpler to use some proper 
approximations so as to be able to carry out a numeric analysis. 

Considering the study [31] , in the NLO case the difference between analytic and stan- 
dard coupling constants can be represented in the form given in (|15p . which is similar to 
the LO one with the additional coefficient equal to 1/2. Note that numerically this NLO 
term is quite analogous to the LO one, since A( ) <C Am (see, for example, [34"]). 

Following the logic expounded in [24J, where it was shown that at the NNLO level the 
effective LO A^ 7 can approximately be taken to be 
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we can apply a simple analytic form in the NNLO as follows H 
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Thus, we propose to use this last expression in the NNLO approximation. 

Since in our analyses Eqs. (0) and were used to derive the strong coupling constant, 
the additional terms in (|14p . (|15p and (|17p can be represented as the functions of the ratio 
A 2 ^/Q 2 , which should be obtained from the corresponding equations 
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5 To have the poles in (fT5j) and (JTTJ) exactly cancelled by those in the perturbative expansions of QCD coupling 
(see, for example, [2]), we keep Q 2 ~ A?^ (i — 0, 1,2) in the denominators of the additional terms (i.e. in the 



last terms of (Tl3|) and (jTTJ) ) and, therefore, above LO we will have additional terms coincinding with those in 
the LO case (with the corresponding replacements Ao — > Ai and Ao — > A) multiplied by the additional factor 



+ [b 2 - b -±jl(a s (Q 2 ))=f3 \n 

In a sense, the replacement quoted in (|14p . (|15p and (|17p is a naive way of doing 
"analytization" ; we apply the latter procedure to the coupling constant itself without 
considering its actual function given by Eqs. ([3]), Q and (fTUj) . Nonetheless, this procedure 
has already been successfully applied in [301 [35] for analyzing the DIS structure functions 
at small x values. With this motivation we would like to investigate its effect in the present 
study as well. Here it will be referred to as a procedure of the naive "analytization". 




3.2 Transition to the canonical form 

To accomplish the procedure of "analytization" more accurately 0, it is convenient for the 
moment (Q 2 ) given in Eq. ([3]) to be represented in the following form (see [36J and 
discussions therein) 



M^(Q 2 ) = R NS (f)-[^ b {Q\Ql) 

where the new "moment" (Q 2 ) is set to be 
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The procedure (I2ip has already been used in [221 [37] , where the Grunberg's effective 
method [38] has been incorporated into the analyses of DIS structure functions. 

Now, the "moment" ji^ s {Q 2 ) has the form close to that obtained for the sum rule, 
because it begins with a s {Q 2 ) rather than a constant. Consequently, the form gets closer 
to that in the difference between the QCD sum rule and its Parton Model value. Follow- 
ing [23], the analytical version of Eq. (|2ip can be obtained by replacing Eq. (|22p for the 
following ones: 
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3 We will call this case an ordinary analytic perturbation theory (APT). 



where Am is the "analytized" m-th power of i-loop QCD coupling p3]. Thus, in the 
APT case the procedure features more complicated functions rather than just the pow- 
ers of some coupling constant. In the one- loop case the Euclidean functions of analytic 
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where L( ) = fn(Q 2 /(A( )) 2 ), while beyond the LO approximation [31\ I24j the transform 
from standard perturbation theory to the APT is slightly modified to assume the following 
form: 
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where the NLO and NNLO results (for i = 1 and 2) are only some approximations. 

To clear up with the meaning of all these formulas, we would like to note one more 
time that in the expressions for Am (Q 2 ) the lower subscript stands for the power of the 
coupling constant, while the upper one is related with the order of the approximation 
considered. Therefore, if the a s -expansion of some variable starts with the first power, as 
is the case at hand, then the power of the last term of the expansion coincides with the 
order of the approximation, i.e. for the last term upper and lower subscripts coincide. 



3.3 Fractional analytic perturbation theory 

Recently, in a series of papers [39], the analytic continuation ([26]) has been extended to 
the noninteger powers of LO L7^,. In the one-loop case, this so-called fractional APT 
(FAPT) gives: 

J_^J__Lw(e-^) 
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is actually the polylogarithm function. 
It is clearly seen that by virtue of 

2 r , ^ ( d \ N z ( d \ N 1 

= — • Ll^W - — - (, s j — , (30) 

all the equations quoted in (|26j) can be shown to be well reproduced. 



'In 139 it was called the Lerch transcendent function. 



Since the product of T-functions satisfy 
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it is easy to obtain the following representation (see Appendix A) : 
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that holds for 9?e(l — i/) > and Me(u) > (in the case at hand these boil down to 
< v < 1). 

Then, following the previous section let's recast Eq, (|28p in the following fashion: 

Lii-„ (A 2 0) /Q 2 

which reproduces the equations in (|26|) with v = 1, 2 and 3, in order. 

Above LO, using quite the same arguments as in the previous subsection we have 
similarly 

(«?W))%W) = (M^y 

which in turn reproduces a set of equations quoted in (f27|) for v = 1, 2 and 3, respectively. 

Note that the Mellin moments (fTUj) used in this case, contain v = 7^(n)/(2/?o) + i 
with i = in the LO, i = 0,1 — NLO, and i = 0,1,2 — NNLO approximations. 
The argument of the polylogarithm function A 2 ^/Q 2 found in Eqs. (I32j) and (|33j) can 
be expressed through the strong coupling constant as was explained in Subsec. 13.11 (see 
Eqs. (fT8l) - (j20l) ). Therefore, the power v lies within the range < v < 4, because < 
7^(n)/(2/?o) < 2 for 2 < n < 10 used in the analyses. The integral representation given 
in Eq. (131 1) is correct only for < v < 1, hence the need to extend it to higher values of 
v. Omitting details of this latter extension H we have 
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To cover an entire range < v < 4, we should consider N = 0,1,2 and 3. For 
N = 0, Eq. (|31|) with the corresponding replacement v — > 5 can be used. For N > 0, it is 
straightforward to obtain 
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and make use of the following formulas valid for < 5 < 1: 
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3 This is done in Appendix A (see Eqs. (|A"3t and (|A"4|| ) for f = iV + (5, where TV = 1, 2, 3. 



i.e. Eqs. ([36]) . ([37|) and ([38]) can be used within the ranges 1<^<2, 2<^<3 and 
3 < v < 4, respectively. 

Note that for 5 = 1 — 5 we have (5 = 1 — 5, 1 + 5 = 2 — 5 and 2 + 5 = 3 — 5. Hence, Eqs. 
(|5B1) . (|37|) and (f38l) can also be used in the ranges < z/ < 1, 1<^<2 and 2 < ^ < 3, 
respectively, if we replace 5 -4 —5 in their r.h.s. It is a strong cross-check of the formulas 
presented in (1361). (pTjl and flggp . 



4 A fitting procedure 

A numeric procedure of data fitting is described in the previous papers [21 [3]. Here we 
just recall some aspects of the so-called polynomial expansion method. The latter was 
first proposed in [40J and further developed in jH]. In these papers the method was based 
on the Bernstein polynomials and subsequently used to analyze data at NLO [121 [25] and 
NNLO level [43} I44j . The Jacobi polynomials for that purpose were first proposed and 
then subsequently developed in [181 H33 120] and used in [221 1251 1251 1331 157] 

With the QCD expressions for the Mellin moments M^{Q 2 ) (k=pQCD,twist2, ...) 
analytically calculated according to the formula in ([2]), the SF F 2 {x, Q 2 ) is reconstructed 
by using the Jacobi polynomial expansion method: 



Nmax 

F 2 k (x,Q 2 )=x a (l-x) b £ (x)$>5 n) (a,/3)M^ 2 (^), 

n=0 j=0 

where 0% b are the Jacobi polynomials and a, b are the parameters to be fit. A condition 
put on the latter is the requirement of the error minimization while reconstructing the 
structure functions. 

As the twist expansion starts to be applicable only above Q 2 ~ 1 GeV 2 the cut Q 2 > 1 
GeV 2 on the data is applied throughout. 

MINUIT program [45] is used to minimize two variables 
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D = dlni^/dlnlnQ 2 . The quality of the fits is characterized by x 2 /DOF for the struc- 
ture function F 2 . Analysis is also performed for the SF slope D that serves the purpose 
of checking the properties of fits (for more details see [2])- 

We use free normalizations of the data for different experiments. For a reference 
set, the most stable deuterium BCDMS data at the value of the beam initial energy 
Eq = 200 GeV is used. With the other datasets taken to be a reference one the variation 
in the results is still negligible. In the case of the fixed normalization for each and all 
datasets the fits tend to yield a little bit worse x 2 > just as in the previous studies. 



5 Results 

Since there is no gluons in the nonsinglet approximation the analysis is essentially easier 
to conduct, with the cut imposed on the Bjorken variable x > 0.25 where gluon density 
is believed to be negligible. 

Here we conduct separate and combined analyses of SLAC, BCDMS and NMC datasets 
obtained with hydrogen and deuterium targets. The cut on x is imposed in a combination 
with those placed on the y = {Eq — E)/Eq variable (£"o and E are initial and final lepton 
energies, respectively) as follows: 

y > 0.14 for 0.3 < x < 0.4 



y > 0.16 for 0.4 < x < 0.5 

y > 0.23 for 0.5 < x < 0.6 

y > 0.24 for 0.6 < x < 0.7 

y>0.25 for 0.7 < x < 0.8 , 



which are meant to cut out those points with large systematic errors. Thus, upon imposing 
the cuts a complete dataset consists of 327 points in the case of hydrogen target and 288 
- deuterium one. The starting point of the QCD evolution is taken to be Qq = 90 GeV 2 . 
This Qq value is close to the average values of Q 2 spanning the corresponding data. From 
earlier studies [21 [3] it follows that it is enough to limit the number of moments to be 
accounted for as in [I9j[20]: N max = 8. Also note that the cut 0.25 < x < 0.8 is imposed 
throughout. 

To reduce the number of parameters, we perform two groups of fits. The first one is 
dealt within a variable-flavor-number-scheme (VFNS) [2j and the and D2 experimental 
datasets analyzed simultaneously. The results for the second group are obtained within a 
fixed-flavor-number-scheme (FFNS) with an active number of flavors n/ = 4 and the H2 
and Z?2 experimental datasets considered separately. 



5.1 VFNS case 

All the results obtained within our reference VFNS [2| and in the cases of "naive", APT, 
and "frozen" modifications of a s are gathered in a set of tables separately for each order 
of perturbation theory approximation and displayed in Figs. 1-3 separately for all these 
three cases. 



Table 1. LO values of the twist-four term h^(x) (with statistic errors given) ob- 
tained in the analysis of the combined H2 + D% dataset within a VFNS and with various 
modifications of a s 



X 


Naive analyt. a B 


APT-inspired a s 


Frozen a s 


Standard a s 


0.275 


-0.121 ± 0.008 


-0.123 ± 0.008 


-0.204 ± 0.011 


-0.271 ±0.012 


0.35 


-0.055 ± 0.007 


-0.055 ± 0.007 


-0.167 ± 0.017 


-0.257 ±0.017 


0.45 


0.119 ± 0.012 


0.119 ± 0.012 


-0.021 ± 0.031 


-0.144 ±0.030 


0.55 


0.422 ± 0.022 


0.422 ± 0.023 


0.211 ± 0.053 


0.051 ±0.049 


0.65 


0.870 ± 0.060 


0.866 ± 0.059 


0.558 ± 0.095 


0.364 ±0.088 


0.75 


1.322 ± 0.117 


1.336 ± 0.112 


0.917 ± 0.152 


0.709 ±0.138 


X z /DOF 


0.93 


0.93 


0.91 


0.94 


X 2 8lop JDOF 


2.30 


2.35 


2.15 


2.60 




0.1474 


0.1474 


0.1409 


0.1400 



Table 2. NLO values of the twist-four term h±(x) (with statistic errors given) ob- 
tained in the analysis of the combined H2 + D2 dataset within a VFNS and with various 
modifications of ct s 



X 


Naive analyt. ot s 


APT-inspired a s 


Frozen a s 


Standard a s 


0.275 


-0.127 ± 0.009 


-0.129 ± 0.007 


-0.183 ± 0.008 


-0.229 ± 0.010 


0.35 


-0.098 ± 0.007 


-0.024 ± 0.009 


-0.149 ± 0.010 


-0.218 ± 0.016 


0.45 


0.014 ± 0.012 


0.187 ± 0.013 


0.010 ± 0.019 


-0.084 ± 0.030 


0.55 


0.172 ± 0.024 


0.506 ± 0.019 


0.215 ± 0.033 


0.098 ± 0.052 


0.65 


0.339 ± 0.057 


0.910 ± 0.045 


0.476 ± 0.065 


0.356 ± 0.093 


0.75 


0.478 ± 0.107 


1.230 ± 0.090 


0.757 ± 0.108 


0.648 ± 0.145 


X^/DOF 


0.85 


0.84 


0.97 


1.02 


~% sloped DOF 


0.82 


0.78 


0.87 


1.20 


Q S (A/|) 


0.1275 


0.1224 


0.1169 


0.1152 



Table 3. NNLO values of the twist-four term h^(x) (with statistic errors given) 
obtained in the analysis of the combined H2 + D2 dataset within a VFNS and with 
various modifications of a s 



X 


Naive analyt. a s 


APT-inspired a s 


Frozen a s 


Standard a a 


0.275 


-0.171 ± 0.006 


-0.196 ± 0.008 


-0.149 ± 0.006 


-0.173 ± 0.017 


0.35 


-0.160 ± 0.008 


-0.152 ± 0.012 


-0.129 ± 0.013 


-0.094 ± 0.020 


0.45 


-0.044 ± 0.018 


0.030 ± 0.022 


-0.007 ± 0.031 


-0.110 ± 0.015 


0.55 


0.085 ± 0.033 


0.269 ± 0.038 


0.116 ± 0.062 


-0.086 ± 0.033 


0.65 


0.221 ± 0.065 


0.551 ± 0.074 


0.218 ± 0.115 


0.085 ± 0.083 


0.75 


0.304 ± 0.100 


0.782 ± 0.116 


0.258 ± 0.169 


0.158 ± 0.105 


X^/DOF 


0.98 


1.02 


0.97 


0.92 


~% sloped DOF 


1.22 


1.17 


1.02 


1.83 


a s (M|) 


0.1151 


0.1125 


0.1163 


0.1159 



From Tables 1-3 it is seen that in all the cases considered x 2 /DOF shows good agree- 
ment between experimental data and theoretical predictions for the SF Mellin moments. 
In LO and NLO cases, the analytic and "frozen" modifications lead to some additional 
improvement of fits, namely, x 2 /DOF in these cases is less than that in the standard 
case. 

Since the quantity xtiope/ DOF is inherently linked with the pQCD aspects to be 
observed in the data, it in this respect is very informative for it strongly varies from one 
approximation to the other thus indicating if there are any effects incompatible with the 
Q 2 dependence in each x-bin assumed. It is seen that it, much like the x 2 /DOF quantity, 
demonstrates similar tendency: in all the cases considered, the analytic and "frozen" 
modifications lead to improvement of fits. Certain improvement of the quality of fits is 
observed at the NLO level. In this approximation even the standard a s case leads to 
reasonable agreement for the slopes; furthermore, for all the infrared modifications it is 
seen that xliope/ DOF < 1. In the NNLO approximation the numbers for the slope in 
the case with a standard a s are not as good but the infrared modifications (especially the 
"frozen" one) lead to better agreement between theoretical and experimental results for 
the former quantity. It is difficult to pin down a reason for such a deterioration of this 
agreement (in the scheme with a standard a s ) when NNLO corrections are added. Because 
this effect is absent in the FFNS case (see the following subsection), we suppose that the 
exact equation given in ([Jj) for the coupling constant is somehow in inconsistency with 
the NNLO expression for the heavy quark thresholds, which is in fact based on certain 
expansions of the coupling constant at the threshold crossing points (see, for example, [2]). 
We plan to study this fine effect elsewhere. 

The QCD coupling decreases from LO through NNLO, which is in perfect agreement 
with other studies (see [36] and references therein); it is seen that the frozen modification 
gives the results closest to those for the standard a s . In the analytic cases the values 
of a s (M|) are higher in the first two orders of perturbation theory, with the maxima in 
the difference Aa s (M§) being 0.0123 and 0.0072 for the naive analytic and APT cases, 
respectively, observed at NLO. Thus, the differences are substantially greater than the 
values of the total experimental error A total a s (M%) = 0.0022, which was obtained in [2] 
by combining statistical and systematic errors in quadrature. 

At the NNLO the APT procedure leads to lower a s (jVf|), though in the naive case 
the a s (M§) value is closer to those obtained for the frozen and standard versions. Nev- 
ertheless, all the NNLO a s (M§) values, except for the APT case, are in good agreement 
within statistical errors, which were found to be /\ stat a s {M 2 z ) = 0.0007 m our previous 
studies [2]. The APT-inspired QCD coupling constant is compatible with the rest within 
the total experimental error. 

From Tables 1-3 and, particularly, Figs. 1-3 it is seen that in the cases of analytic and 
"frozen" coupling constants the twist-four corrections are larger compared to those in the 
case of a standard perturbative coupling constant, thus confirming the results obtained 
in [23]. For example, at x ~ 0.75 the higher-twist corrections (HTCs) in the standard 
case are about twice as less that those obtained in the cases of analytic modifications in 
all orders considered. In the "frozen" case, the HTC values are compatible with those in 



the standard case within statistical errors. □ 

The difference in the values of HT parameter for the naive analytic and frozen variants 
becomes moderate at the NLO level unlike the APT case. For the infrared modifications 
the HT terms are large and in the frozen and APT cases are compatible with the LO ones. 

In the NNLO approximation the situation changes drastically. HT terms for all the 
cases, except for the APT case, are comparable with each other. HTCs in the latter case 
are still higher but they also strongly suppressed to compare with them in LO and NLO 
cases. Thus, at the NNLO level the twist-four corrections appear to be small for all the 
cases considered above. 

The values for parameters in the parameterizations of the parton distributions (see 
Eq. @) for the cases corresponding to different coupling constant modifications are given 
in Table 4. From this table it is seen that the parameters in all the cases considered are 
close to each other. 

Moreover, the parameter b, which is responsible for the large x fall of the PDF shape, 
is found to be around 4, that is in agreement with the quark counting rules |48j . if we 
take into account the rise of b at large Qq values used here (see Eq. (|B2p ). 



Table 4. The values of PDF parameters obtained within a VFNS and with various 
modifications a s in the case of the combined H2 + D2 dataset analyzed 





Naive analyt. a a 


APT-inspired ct s 


Frozen a s 


Standarc 


a s 


Par. 


LO 


NLO 


NNLO 


LO 


NLO 


NNLO 


LO 


NLO 


NNLO 


LO 


NLO 


NNLO 


A(H 2 ) 


2.43 


2.65 


2.54 


2.46 


2.54 


2.51 


2.48 


2.63 


2.55 


2.53 


2.62 


2.57 


b(H 2 ) 


3.73 


4.06 


4.16 


3.72 


4.07 


4.14 


3.71 


4.03 


4.16 


3.71 


4.02 


4.13 


d(H 2 ) 


5.26 


4.72 


6.13 


5.13 


5.61 


6.29 


5.06 


4.95 


5.98 


4.92 


5.07 


5.60 


A(D 2 ) 


2.29 


2.64 


2.43 


2.32 


2.52 


2.45 


2.31 


2.57 


2.48 


2.32 


2.52 


2.50 


b(-Da) 


3.77 


4.09 


4.21 


3.76 


4.11 


4.19 


3.75 


4.07 


4.20 


3.75 


4.07 


4.15 


d(D 2 ) 


3.07 


2.79 


3.72 


2.97 


3.21 


3.66 


2.99 


2.96 


3.57 


2.95 


3.11 


3.39 



5.2 FFNS case 

For comparison let's present the values of PDF parameters obtained within the FFNS 
(nf = 4) in the analyses of the hydrogen data for the versions of the strong coupling 
constant discussed above, with the addition for the case of FAPT-inspired modification 
of a s . 

Here we have no threshold transitions for QCD coupling and PDF Mellin moments 
and, hence, are able to consider the hydrogen and deuterium datasets separately. From 
Table 5, we see that all the parameters in the parameterizations of the parton distributions 
(see Eq. ([3D) for the cases corresponding to different coupling constant modifications are 
close to each other in all the cases considered. As was already discussed in the previous 
subsection, the parameter b ~ 4, that is in agreement with the quark counting rules [48J , 
if we take into account its certain rise at large Qq values. 



Table 5. Tiie values of PDF parameters obtained within the FFNS (nf = 4) for 
various modifications of a s in the case of H2 and D2 datasets analyzed separately 





Naive analyt. a s 


APT a s 


Frozen a s 


FAPT a s 


Par. 


LO 


NLO 


NNLO 


LO 


NLO 


NNLO 


LO 


NLO 


NNLO 


LO 


NLO 


NNLO 


A(H 2 ) 


2.30 


2.58 


2.49 


2.30 


2.52 


2.42 


2.32 


2.57 


2.49 


1.68 


2.55 


2.47 


b(H 2 ) 


3.74 


4.06 


4.18 


3.74 


4.08 


4.19 


3.74 


4.04 


4.18 


3.74 


4.05 


4.18 


d(H 2 ) 


5.63 


5.30 


6.46 


5.63 


5.82 


6.91 


5.70 


5.36 


6.38 


5.90 


5.50 


6.54 


A(D 2 ) 


2.42 


2.75 


2.60 


2.42 


2.57 


2.51 


2.38 


2.66 


2.62 


2.34 


2.65 


2.58 


b(D 2 ) 


3.73 


4.05 


4.19 


3.73 


4.08 


4.20 


3.72 


4.02 


4.17 


3.72 


4.03 


4.17 


d(D 2 ) 


2.61 


2.47 


3.15 


2.61 


2.98 


3.42 


2.72 


2.61 


3.10 


2.82 


2.64 


3.18 



9 For the results to be discernible, Figs. 1-4 and 6-7 contain statistical errors for the infrared modifications 
only. The magnitude of the errors is on par in the standard case and is not shown. 



By using the case of "naively analytized" a s it is demonstrated in Fig. 5 that the 
values of HTCs obtained for the n/ = 4 and n/ = 5 cases in the FFNS are very much 
alike. 

In order to be able to assess the difference among the cases with frozen, FAPT and 
standard versions for the strong coupling constant let's present the tables with the HT 
values obtained within FFNS (n/ = 4). 

Table 6. LO values of the twist-four term obtained within the FFNS (rif = 4) and 
with various modihcations of a s in the separate analyses of H2 and D2 datasets 





FAPT-inspircd a a (± stat) 


Frozen », 


(± stat) 


Standard a s (± stat) 


X 


h 4 (x) for H2 


h 4 (x) for D2 


h 4 (x) for H2 


h 4 (x) for £>2 


h 4 (x) for H2 


h 4 (x) for D2 


0.275 


-0.221±0.011 


-0.206±0.010 


-0.183 ± 0.012 


-0.167 ± 0.015 


-0.235±0.012 


-0.229±0.011 


0.35 


-0.187±0.014 


-0.134±0.009 


-0.160 ± 0.018 


-0.104 ± 0.022 


-0.232±0.021 


-0.192±0.010 


0.45 


0.002±0.023 


0.072±0.016 


-0.023 ± 0.037 


0.055 ± 0.042 


-0.130±0.038 


-0.072±0.018 


0.55 


0.332±0.034 


0.462±0.026 


0.189 ± 0.065 


0.330 ± 0.074 


0.049±0.065 


0.158±0.029 


0.65 


1.001±0.063 


1.078±0.068 


0.610 ± 0.117 


0.701 ± 0.132 


0.455±0.111 


0.495±0.067 


0.75 


2.031±0.131 


1.701±0.142 


1.177 ± 0.207 


0.913 ± 0.206 


1.003±0.182 


0.677±0.128 


X z /DOF 


0.98 


0.77 


0.94 


0.75 


0.98 


0.78 


^ sloped 


1.47 


2.23 


1.25 


1.93 


1.58 


2.45 


a s (M|) 


0.1394 


0.1426 


0.1387 


0.1388 


0.1376 


0.1383 



Table 7. NLO values of the twist-four term obtained within the FFNS (rif = 4) and 
with various modifications of a s in the separate analyses of H2 and D2 datasets 





FAPT-inspircd a a (± stat) 


Frozen a t 


(± stat) 


Standard a s (± stat) 


X 


h A (x) for H 2 


h A (x) for D 2 


h 4 (x) for H 2 


h 4 (x) for D 2 


h 4 (x) for H 2 


h 4 (x) for D 2 


0.275 


-0.205±0.012 


-0.209±0.019 


-0.154 ± 0.011 


-0.147 ± 0.014 


-0.200±0.012 


-0.206±0.015 


0.35 


-0.194±0.017 


-0.165±0.026 


-0.147 ± 0.016 


-0.099 ± 0.021 


-0.219±0.019 


-0.189±0.024 


0.45 


-0.067±0.030 


-0.027±0.039 


-0.058 ± 0.037 


0.013 ± 0.044 


-0.169±0.041 


-0.122±0.048 


0.55 


0.160±0.044 


0.252±0.052 


0.069 ± 0.067 


0.204 ± 0.078 


-0.078±0.072 


0.018±0.083 


0.65 


0.635±0.071 


0.702±0.083 


0.317 ± 0.118 


0.432 ± 0.139 


0.153±0.123 


0.207±0.143 


0.75 


1.512±0.130 


1.243±0.137 


0.723 ± 0.204 


0.509 ± 0.210 


0.534±0.205 


0.251±0.212 


x '/DOF 


0.91 


0.70 


0.89 


0.71 


0.92 


0.73 


X 2 slop JDOF 


1.08 


1.35 


0.92 


1.27 


1.23 


1.68 




0.1220 


0.1261 


0.1200 


0.1201 


0.1192 


0.1199 



Table 8. NNLO values of the twist-four term obtained within the FFNS (n f = 4) 
and with various modihcations of a s in the separate analyses of H2 and D2 datasets 





FAPT-inspircd a a (± stat) 


Frozen », 


(± stat) 


Standard a s (± stat) 


X 


h 4 (x) for H2 


h 4 (x) for D2 


h 4 (x) for H2 


h 4 (x) for D2 


h 4 (x) for H2 


h 4 (x) for D2 


0.275 


-0.168±0.010 


-0.175±0.013 


-0.119 ± 0.011 


-0.116 ± 0.012 


-0.158 ± 0.020 


-0.161 ± 0.012 


0.35 


-0.181±0.015 


-0.155±0.020 


-0.112 ± 0.010 


-0.072 ± 0.017 


-0.166 ± 0.021 


-0.146 ± 0.017 


0.45 


-0.130±0.037 


-0.095±0.045 


-0.049 ± 0.018 


0.012 ± 0.042 


-0.156 ± 0.036 


-0.108 ± 0.043 


0.55 


-0.056±0.065 


0.023±0.079 


0.007 ± 0.031 


0.129 ± 0.080 


-0.157 ± 0.079 


-0.035 ± 0.080 


0.65 


0.126±0.111 


0.162±0.133 


0.106 ± 0.058 


0.209 ± 0.145 


-0.061 ± 0.127 


0.022 ± 0.142 


0.75 


0.419±0.178 


0.163±0.191 


0.240 ± 0.116 


0.088 ± 0.213 


0.049 ± 0.209 


-0.117 ± 0.207 


X z /DOF 


0.89 


0.70 


0.87 


0.68 


0.89 


0.71 


X 2 slop JDOF 


0.97 


1.22 


0.68 


0.85 


0.97 


1.28 


a s (M 2 z ) 


0.1183 


0.1197 


0.1180 


0.1184 


0.1176 


0.1176 



Just like in the previous subsection it is seen from Tables 6-8 that in all the cases at 
hand, we have good agreement between the experimental data and theoretical predictions. 
Once again, in all the cases considered, excluding LO FAPT, the analytic and frozen 
modifications lead to slight improvement of fits, that is y^/DOF and )^ slope /DOF are 
found to be smaller in these cases than those for the standard a s , and the QCD coupling 
constant decreases when we move from LO through NNLO. 

Also note that contrary to what was observed in the previous subsection here the 
quantity xliope/ DOF steadily decreases when we proceed step by step from LO to NLO 



and then to NNLO level. Respectively, here the smaller values of xti op e/ DOF are observed 
for the "frozen" case as well. 

The values of coupling constants are very similar, especially for the "frozen" and 
standard versions. In the analytic cases, the central values of a s (M|) are little higher, 
mostly in the first two orders of perturbation theory. 

In the NNLO, all the a s (M§) values, except for the analytic one, are in good agreement 
within statistical errors. As earlier, the analytic QCD coupling is in agreement with the 
rest within the total experimental error. 

Similarly to the previous subsection we note in Figs. 4, 6, 8 appreciable difference 
between HTCs obtained in the standard and analytic cases in the first two orders of 
perturbation theory. The difference in the cases of a "frozen" and standard a s is not as 
large; the values of HT terms are in agreement within statistical errors. 

In the NNLO, the situation is over again changed considerably. In all the cases of in- 
frared modifications considered, the HTCs are small. Similarly to the previous subsection 
they are not compatible with zero at x ~ 0.75, while the HT terms in standard QCD are 
compatible with zero and, at the same time, in agreement with all the cases considered 
within statistical errors. 

Then, at NNLO we can see agreement between standard QCD and its infrared modi- 
fications for QCD coupling a s (M|), as well as for the respective HTCs (as a rule) within 
statistical errors. 

6 Conclusions 

The pattern of separating perturbative QCD and HT corrections may be different in 
different orders of perturbation theory, as well as in some resummations based on several 
first orders, and for certain modifications of the strong coupling constant as well. In 
the present paper, we studied consequences of the infrared modifications of the QCD 
coupling constant — the so-called "frozen" and "analytized" versions. In the last case we 
considered three different options: 

• a simple modification of the strong coupling constant [H] without rearrangement of 
a perturbation series; 

• an application of the ordinary analytic perturbation theory (see |24j ) to the "mo- 
ments" fJ-n S (Q 2 ) given in (1221) ; 

• impact of the fractional analytic perturbation theory [39] applied directly to the 
Mellin moments M^ S (Q 2 ). 

To test all these modifications, the Jacobi polynomial expansion method developed in \18\ 
fl~9j [20] was used to perform analysis of (^-evolution of the DIS structure function F2 
by fitting all existing to date reliable fixed-target experimental data that satisfy the cut 
x > 0.25. To the best of our knowledge, this is the first application of the FAPT results 
to the fits of the DIS structure functions. 

From the results obtained we conclude that in the first two orders of perturbation 
theory the value of the QCD coupling constant at the reference point a s (M|) is larger with 
respect to that in the standard case, for all versions of the analytic perturbation theory. 
An increase in the central values of a s (M§) is smaller in the NNLO approximation, much 
like the case of the "frozen" version of a s . Nevertheless, in all the cases the a s (M^) values 
agree mostly within total experimental errors. For the NNLO case, the results are as a 
rule compatible between each other within statistical errors. Also, we note that within 
statistical errors only (in our case A stat a s (M^) = 0.0007 [2]) all the NNLO results for 



10 Since the magnitude of HTCs is a bit higher for the hydrogen dataset (see, for example, our recent paper 
we consider only this case. 



a s (M§) in the FFNS case (see Table 8) are too compatible with the world average value 
for the coupling constant presented in the review 

a s {M 2 z ) = 0.1184 ±0.0007. 

It is also observed that the values of the twist-four corrections in the first two orders 
of perturbation theory are higher than those in the standard case, particularly for all 
versions of the analytic perturbation theory. This observation is in complete agreement 
with recent studies [24J of the Bjorken sum rules, where it was shown that there is a 
reduction of higher HT terms, starting with the twist-six ones. Unfortunately, here we 
are not able to study the twist-four and twist-six corrections simultaneously since their 
contributions are strongly correlated. 

In the NNLO, all the cases of infrared modifications of the QCD coupling feature 
nonzero although rather small twist-four corrections. Whereas in the case of the standard 
QCD approach, the HT terms are close to zero at large x values. However, the NNLO 
HTCs for all the cases considered are compatible between each other within statistical 
errors. In principle, the main difference between standard QCD and its analytic and 
"frozen" modifications is in the strong suppression of the higher twist corrections in NLO 
and NNLO orders of perturbation theory, respectively. 

What is interesting to look for further in the study is the consideration of the com- 
bined nonsinglet and singlet analyses using the DIS experimental data within an entire 
x region, as well as an application of certain resummation-like Grunberg effective charge 
methods [38] (as was done in |37| in the NLO approximation) and the "frozen" [13] 1^1 and 
analytic [3TJ versions of the strong coupling constant (see [351 GEH 121] for recent studies 
in this direction). The effect of N 3 LO corrections would also seem to be important to 
account for in the subsequent investigations, as well as an extension of the FAPT model 
results for a* to the VFNS case. 
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8 Appendix A 

Here we present a simple evaluation of the equation given in (I3ip . 
Using Sec. II. 5 in [46J, we have 

oo 

Ul -W = fff^J^ (»<!)■ (AD 

Upon substituting t = ln(l/£), we obtain for the r.h.s. of Eq. 

l 

I'M 







(A2) 



11 It should be mentioned that this analysis was carried out over the data coming from the various experiments 
and in different orders of perturbation theory, i.e. from NLO up to N 3 LO. 

12 There are a lot of "frozen" versions of the strong coupling constant (see, for example, the list of references 
in [30]). 



In order to extend Eq. (|A2p to the values of v > 1, it is useful to consider a series 
representation (I29p for u = N + 5: 



E 



N+S-l m 



m Z 
m=l x "~/ 171=1 



j x N oo 

dz I ^ 



N 



z^- z ) Li ,_,(•-). (A3) 



Taking the representation (|A2|) for the r.h.s. of (jA3j), we have 

i 



i'V 



r(i-j) ji-^ U 



(0 < 5 < 1) . (A4) 



9 Appendix B 

Consider the large x asymptotic [47J of SF F2(x,Q 2 ) (for simplicity we restrict ourselves 
to the LO approximation) 



.2, . , , ^b(a{ 0) (Q 2 )) 



F 2 (x,Q')~(l-x) b 



(Bl) 



where 



bo-dHa^iQ 2 )), d 



16 

3ft ' 



and bo is some constant which can be obtained from the quark counting rules 
The Q 2 -dependent part in the r.h.s. of (|Blj) can be represented as follows: 



(l_ x )-dln(4 0) (Q 2 )) 



-dln(l-z) 



Performing the "analytization" procedure and using (|14p . we have 

N -dln(l-a;) 



A ?o) 



(0), 



<4 0) (Q 2 ) 



-dln(l-x) 



1 + 



dln(l-x) A ? ( 



(0) 



(B2) 

From this expression it follows that the power corrections, which were generated by the 
"analytization" of the coupling constant, have the opposite sign as compared to the twist- 
four corrections at large x and, moreover, demonstrate a different asymptotic behavior. 
These corrections (taken with the additional sign "— ") increase like ln(l — x) while the 
twist-four corrections behave like 1/(1 — x) (see [27] and discussions therein). This dif- 
ference in sign between these corrections leads to the larger twist-four corrections in the 
present analysis than those given in [2]. 
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Figure 1: Comparison of the HTC parameter h±{x) obtained in LO, NLO and NNLO for 
hydrogen data (the bars show statistical errors) between our reference VFNS with a standard 
perturbative a s [2] and that with naive analytic a s . 
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Figure 2: Comparison of the HTC parameter h±(x) obtained in LO, NLO and NNLO for 
hydrogen data between a VFNS and that with APT-inspired a s . 
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Figure 3: Comparison of the HTC parameter h^(x) obtained at LO, NLO and NNLO for 
hydrogen data within a VFNS between the cases with a standard perturbative and "frozen" 
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Figure 4: Comparison of the HTC parameter h^x) obtained in LO, NLO and NNLO for 
hydrogen data within a FFNS (n/ = 4) between the cases with a standard perturbative and 
naively analytized ot s . 
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Figure 5: Comparison of the HTC parameter h^{x) obtained in LO, NLO and NNLO for 
hydrogen data within a FFNS between the cases with nf — 4 and rif = 5 and the naive 
analytic a s . 
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Figure 6: Comparison of the HTC parameter h^{x) obtained in LO, NLO and NNLO for 
hydrogen data within a FFNS (n/ = 4) between the cases with a standard perturbative and 
APT- inspired a s . 
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Figure 7: Comparison of the HTC parameter h±{x) obtained at LO, NLO and NNLO for 
hydrogen data within a FFNS (n/ = 4) between the cases with a standard perturbative and 
"frozen" a,. 
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Figure 8: Comparison of the HTC parameter h^{x) obtained at LO, NLO and NNLO for 
hydrogen data within a FFNS (n/ = 4) between the cases with a standard perturbative and 
FAPT-inspired ot s . 



